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بهینه‌سازی الگوی کشت در راستای مدیریت منابع آب‌های زیرزمینی با استفاده از الگوریتم‌های 
5 و ۳50 


مهدی کماسی "*- امیر علیزاده فرد۲ 
تاریخ دریافت: ۱۳۹۹/۱۲/۰۵ 


تاریخ پذیرش: ۱۴۰۰/۰۲/۱۴ 


چکیده 


پایش‌های میدانی بیانگر کاهش تراز آبخوان‌های دشست سیلاخور در سال‌های اخیر هستند که این موضوع لزوم بهره‌برداری بهینه از منابع محدود 
آب‌های زیرزمینی را می‌رساند. در همین راستا باتوجه‌به سهم عمده کشاورزی در مصرف آب بهینه‌سازی الگوی کشت از مهم‌ترین عوامل مدیریت منابع 
آب و مقابله با کم‌آبی محسوب می‌شود. در این پژوهش عملکرد توابع جریمه مختلف به‌منظور تعیین الگوی کشت بهینه در دو الگوریتم تکامل تفاضلی و 
ازدحام ذرات بررسی شده است. در این راستا در گام نخست بارندگی سال‌های زراعی ۱۳۹۹-۱۴۰۰ ۱۳۰۰-۱۳۰۱ و ۱۴۰۱-۱۴۰۲ در دشت سیلاخور با 
استفاده از برنامه‌ریزی ژنتیکی و مدل میانگین متحرک خود هم‌بسته یکپارچه فصلی پیش‌بینی شده و باتوجه‌به تغذیه ناشی از بارندگی و میزان برداشت 
در سال‌های گذشته. برای هرسال زراعی ۱۰۰ سناریو برداشت تعیین شد. سپس الگوی کشت بهینه متناسب با هر سناریو برداشت با بهره‌گیری از 
برنامه‌ریزی خطی به دست آمد. حداکنرسازی درآمد کشاورزان به‌عنوان تابع هدف» و آب و زمین در دسترس به‌عنوان محدودیت‌های برنامه‌ریزی خطی در 
نظر گرفته شدند. در گام دوم عملکرد توابع جریمه مختلف در دو الگوربتم تکامل تفاضلی و ازدحام ذرات برای رسیدن به پاسخ‌های حاصل از برنامه‌ریزی 
خطی مورد بررسی قرار گرفت. طبق نتایج این بررسی به طور میانگین الگوریتم ازدحام ذرات در ۸۳/۶ سناریوها سریع‌تر به الگوی کشت بهینه می‌رسد. 
همچنین اعمال تابع جریمه پویای طبقه‌بندی شده تعداد تکرارهای الگوریتم تکامل تفاضلی را ۱۱/۶ و الگوریتم ازدحام ذرات را /۳ نسبت به استفاده از 
جریمه ساکن کاهش می‌دهد. در مجموع می‌توان دریافت که با میزان برداشت فعلی» در صورت پیروی از الگوی کشت بهینه پیشنهاد شده برای چهار 
محصول عمده زراعی دشت سیلاخوره یعنی گندم» جو شلتوک و ذرت علوفه‌ای. درآمد کشاورزان منطقه ۱۸/۶ افزایش می‌بابد. همچنین به کارگیری تابع 


جریمه پویای طبقه‌بندی شده در الگوریتم ازدحام ذرات به‌منظور بهینه‌سازی الگوی کشت پیشنهاد می‌شود. 


واژه‌های کلیدی: آب‌های زیرزمینی» الگوریتم‌های فراکاوشی برنامه‌ریزی خطی بهینه‌سازی, دشت سیلاخور 


مقدمه 

وقوع خشکسالی‌های پی درپی در کنار افزایش نیازهای آبی و عدم 
مدیریت صحیح منابع آب» باعث ایجاد بحران آبی شده است که 
پیامدهای مختلف محیطزیستی و اقتصادی را به‌دنبال دارد. از آن جا 
که این بحران بر بخش‌های مختلف جامعه تاثیر می‌گذارد. بنابراین 
بررسی راه‌های مقابله با آن نیزه امری لازم و ضروری است. در زمینه 
شناسایی علل خشکسالی رویکرد واحد و روشنی وجود ندارد. برخی 
تغییر اقلیم گسترده را به‌عنوان عامل اصلی بحران کم‌آبی؛ در ایران 
معرفی می‌کنند (۲۳ و ۳۰) و در مقابل برخی دیگر علت آن را برداشت 


۱ و ۲- به‌ترتیب دانشیار و دانشجوی کارشناسی ارشد گروه مهندسی عمران» 
دانشگاه آیت اه العظمی بروجردی» بروجرد» ایران 


بی‌روبه انسان از منابع آبی و مدیریت نادرست منابع آب می‌دانند (۲ و 
٩‏ نتایج پژوهش‌ها حاکی از آن است که در دشت سیلاخور عوامل 
انسانی نسبت به تغیبرات اقلیمی تاثیر بیشتری در کاهش تراز آب‌های 
زیرزمینی دارند (۱۲). علاوه بر خشسکس‌الی. تغییر الگوی کشت به 
سمت محصولات آب‌بر نیز بحران آب را در سال‌های اخیر به نخستین 
پدیده بحرانی در سطح جامعه تبدیل کرده‌است که این بحران بر روی 
بخش کشاورزی به‌عنوان بزرگ‌ترین مصرف کننده آب تاثیر مستقیم 
می‌گذارد. آب از مهم‌ترین اجزاء برای تولید محصولات کشاورزی است 
و مدیریت منابع ارزشمند و محدود آن از زمینه‌های گسترش کشاورزی 
می‌باشد. تولید محصولات کشاورزی به طور متوسط ۷۰/۶ از کل 
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(#- نویسنده مسئول: 


۳۵۰ تشریه آب و خاک. جلد ۰۲۵ شماره ۳. مرداد - شهریور ۱۴۰۰ 


| 
راهکارهای مناسب جهت کاهش هدررفت آب و دستیابی به توسعه 
پایدار بهینه‌سازی تخصیص است (۱۸). در این راستا باید به اهداف 
اقتصادی که شامل درآمد کشاورزان می‌شود نیز توجه نمود. به‌طورکلی 
هدف محیطزیستی برداشت کمتر از آب‌های زیرزمینی و هدف 
اقتصادی افزایش درآمد کشاورزان با یکدیگر در ناسازگاری هستند. در 
بررسی مسئله برنامه‌ریزی مصرف آب در بخش کشاورزی باید تعادل 
مناسبی میان این اهداف ایجاد نمود که در این زمینه رویکردهای 
متفاوتی ارائه گردیده است به‌عنوان نمونه کایر و همکاران (۱۰) با 
استفاده از الگوی برنامه‌ریزی خطی" با هدف حداکثر کردن درآمد 
خالص و صرفه‌جویی در مصرف آب. در پنجاب, الگوی کشت بهینه را 
تعیین کردند. آن‌ها نشان دادند که با اضافه کردن پنبه به الگوی کشت 
منطقه در مصرف آب ۲۶/۵۵/۰ صرفه‌جویی شده و درآمد /۴ افزایش 
اه رای ایا هی رز وتات 
فازی به بررسی مدیریت منابع آب پرداختند. نتایج آن‌ها حاکی از آن 
است که مجموعه‌های فازی قادرند عدم قطعیت دو جانبه‌ای در 
مدل‌ها به‌وجود آورند و تصمیم‌های ارزشمندی مطابق با اولویت‌ها و 
شرایط عملی تصمیم گیرندگان با استفاده از روش حل به‌دست می‌آید. 
به‌منظور بهینه‌سازی الگوی کشت در مشهد. مجیدی و همکاران (۱۵) 
برنامه‌ریزی خطی را با هدف کاهش مصرف آب بکار بردند. آن‌ها 
محدودیت‌های مختلفی مانند مسطح زیرکشسته نیروی کاره سرماید 
بازده برنامه‌ای» ماشین آلات کود شیمیایی. سم شیمیایی و تناوب 
زراعی را در مدل‌سازی خود درنظر گرفتند. براس‌اس نتایج اين 
پژوهشگران. محصولاتی با مصرف بالای آب و بازده برنامه‌ای کمتر, 
در الکو کشت پویته خزقاض فقوت سبالقادر و هبکاران (۱) برا 
بهینه‌سازی الگوی کشت از برنامه‌ریزی ریاضی با هدف‌های کاهش 
مصرف آب. حداکثرسازی بازده سالانه و تخصیص کارآمد زمین‌های 
کشاورزی استفاده کردند. آن‌ها میزان واردات و صادرات محصولات 
کشاورزی را نیز در مدل‌سازی خود درنظر گرفته و با مقایسه دو الکوی 
کشت فعلی و بهینه نان دادند در الگوی کشت بهینه. مصرف آب 
کاهش و سود خالص برآورد شده در حدود ۲/۴۲ میلیارد دلار در 
سال افزایش می‌یابد. در الگوی کشت بهینه به محصولات غلات» 
وه هو یپ ۱ ۵ ین 
اختصاص می‌یابد. بیرهانو و همکاران (۳) نیز با استفاده از برنامه‌ریزی 
خطی الگوی کشت بهینه را در کشسور اتیوپی تعیین کردند. در اين 
پژوهش آن‌ها محدودیت آب و زمین را در چهار سناریو مختلفه و 
هدف‌ها را حداکثرسازی بازده وحداکثرسازی تولید تعریف کردند. نتایج 
اين مطالعه نشان داد بهره‌وری آبباری ۰۴۸/۶ بهره‌وری آب فیزیکی 


(ّ) عمتصصصهه۳۳0 تهمصنا -1 
صممتعصن۲ نالفوع۳ -2 


۲ تا ۱۶۸ بهره‌وری آب اقتصادی ۲۰۵ تا ۲۴۱ عملکرد ۱۰۸ تا 
و سود متوسط ۱۵۲ تا ۲۰۸ از طریق بهینه‌سازی الگوی 
کشت ارزش اقتصادی بالاتری خواهد داشت. در پژوهشی دیگر در 
مصر اسامه و همکاران (۲۰) با درنظر گرفتن میزان استراتژیک بودن 
هر محصول و استفاده از برنامه‌ریزی خطی, الگوی کشت بهینه را 
تعبین کرده و با الگوی کشت فعلی مقایسه کردند. آن‌ها به این نتیجه 
تسیا که در اگوی کفنت یه سح کشت تم خی ماند بیان 
سیر نخود و عدس به‌دلیل استراتژیک نبودن آن‌هاء کاهش می‌بابد. 
درمقابل سطح کشت محصولاتی مانند گندم؛ برنج» ذرت و محصولات 
قندی که ارزش غذایی بالایی دارند و به‌عنوان محصولات استراتژیک 
تعریف می‌شوند» تقریبا ثابت باقی می‌ماند. همچنین مناطق اختصاص 
يافته به محصولاتی که بازده خالص بالایی دارن مانند گوجه فرنگی 
بطور قابل توجهی افزایش خواهد یافت. 

رویکرد دوم در زمینه تخصیص بهینه و تعیین الگوی کشت 
بهره‌گیری از الگوریتم‌های فراکاوشی است؛ به‌عنوان نمونه راجو و 
کومار (۲۴) در راستای افزايش سود یک پروژه آبیاری. عملکرد 
الگوریتم ژنتیک را در برنامه‌ریزی آبیاری و توسعه الگوی کشت بررسی 
کردند. در این مطالعه مجموع حداکثر عملکرد نسبی کل گیاهان در 
مناطق تحت آبیاری به‌عنوان تابع هدف در نظر گرفته سد. 
محدودیت‌های مدل نیز شامل معادله پیوستگی, نیاز آب تنوع 
محصول و محدودیت ذخیره بود. همچنین برای انحراف مسئله از 
خالت. منطقی به غیرمنطقی از تایح جریمه" اسستفاده ند آن‌ها تشسان 
دادند نتایج بهینه به‌دست آمده از الگوریتم ژنتیک و برنامه‌ریزی خطی 
یکسان می‌باشد. در پژوهشی دیگر به‌منظور بررسی شرایط مختلف 
بهره‌برداری از مخازن چندمنظوره در هندوستان ردی و کومار (۲۵) 
الگوریتم ازدحام ذرات مبتنی بر نخبه گرایی" را بکار بردند. اين 
پژوهشگران دخیره‌سازی اولیه. تخلیه» سرریز و تبخیر را در سمت 
مخزن, و تخصیص آب تبخیر و تعرق واقعی و رطوبت اولیه خاک را 
در سمت مزرعه» به عنوان خروجی‌های مدل در هر دوره ده روزه و 
برای هر محصول تعریف کردند. مدل آن‌ها تصميم‌گیری برای 
تخصیص بهینه آب را تسهیل می‌کند. یوسف دوست و همکاران (۳۳) 
از الگوریتم ژنتیک به‌منظور تعیین سطح زیر کشت بهینه برخی از 
محصولات کشاورزی در قزوین» تحت شرایط متفاوت آب و هوایی 
استفاده کردند. بر اساس نتایج این پژوهش» سود حاصل از الگوی 
کشت بهینه در شرایط آب و هوایی نرمال. مرطوب. خشک و گرم و 
خشک افزایش چشمگیری خواهد داشت. از سویی دیگر سیاسر و هنر 
(۲۸) نیز به بهینه‌سازی تخصیص آب و الگوی کشت محصولات غالب 
کشاورزی با بکارگیری الگوریتم ژنتیک در دشت سیستان پرداختند. 
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کماسی و همکاران. بهینه‌سازی الگوی کشت در راستای مدیریت منابع آب‌های زیرزمینی با استفاده از الگور بتم‌های ۱۳ و ۳۹0 ۳۵۱ 


مدل آن‌ها علاوه‌بر الگوی کشت بهینه. سطح زیر کشت محصولات را 
نیز با لحاظ کم‌آبیاری‌های مجاز در کل دوره رشد. به‌صورت همزمان» 
تعیین می‌کند. همچنین اثرمحدودیت آب در دوره‌های مختلف رشد بر 
سطح زیر کشت محصولات مختلف در این مطالعه بررسی شده است. 
تتایج این بررسی حاکی از آن اسست که گندم و هندوانه به‌ترتیب 
محصولات اقتصادی کشت اول و دوم زراعی هستند. 

استفاده از برنامه‌ریزی خطی برای تعیین الگوی کشت در حقیقت 
حل یک مسئله بهینه‌سازی مقید" خطی است. برای حل مسائل مقید 
با اسستفاده از الگوریتم‌های بهینه‌سازی نیز روش‌های مختلفی وجود 
دارد. در این پژوهش پس از پیش‌بینی میزان بارندگی سه سال آینده و 
محاسبه تغذیه ناشی از آن» محدودیت‌های برداشت در ۱۰۰ سناریو 
مختلف برای هرسال تعریف شده و دو رویکرد استفاده از برنامه‌ریزی 
خطی و استفاده از الگوریتم‌های بهینه‌سازی برای تعیین الگوی کشت 
بهینه محصولات عمده زراعی دشت سیلاخور بررسی و مقایسه 
شده‌اند. همچنین عملکرد سه نوع تابع جریمه ساکن " پوبا و پوبای 
طبقه‌بندی " شده در دو الگوریتم بهینه‌سازی تکامل تفاضلی" و ازدحام 


ذرات " مورد بررسی قرار گرفته است. 


منطقه موردمطالعه 

محدوده مطالعاتی دشت سیلاخور پیرامون شهرهای بروجرد و 
دورود بزرگ‌ترین زمین هموار استان لرستان بوده که با مساحت 
۰۱ کیلومترمربع» در غرب ایران و شمال شرقی‌ترین بخش این 
استان واقع شده است. ۸۱٩‏ کیلومتر مربع این محدوده را دشست و 
۲ کیلومتر مربع آن را ارتفاعات تشکیل داده‌اند. دشت سیلاخور در 
شمالی‌ترین قسمت حوزه آبریز کارون بزرگ واقع شده و از اقلیم سرد 
و بارش مناسب برخوردار است؛ اما مطالعات میدانی در سال‌های اخیر 
کاهش تراز آب‌های زیرزمینی این منطقه را نشان می‌دهد. بر اساس 
گزارش‌های شرکت آب منطقه‌ای لرستان, مخازن آب زیرزمینی دشت 
سیلاخور با ۶۲ میلیون متر مکعب کسری مواجه هستند که یکی از 
عوامل آن تغییر الگوی کشت به سمت کشت محصولات آب‌بر 
مر یادف کل ۱ موفیه دفتتو رات ان رشان و ری کون 
نشان می‌دهد. 


مواد و روش‌ها 


در اين پژوهش به‌منظور تعیین الگوی کشت بهینه محصولات 


1- 0089021260 ۸ 

2- ٩02116 ۷ 
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عمده زراعی دشت سیلاخور با هدف حداکثرسازی درآمد کشاورزان در 
سته ال آینده.ه محتودیت‌های آب و زفین در دستنرس» از عیران 
تغذیه ناشی از بارندگی به‌عنوان معیاری برای مش‌خص کردن بازه 
استخراج آب استفاده شد. به‌منظور پیش‌بینی بارندگی از الگوهای 
سری زمانی 5۸/1/۸" و برنامه‌ریزی ژنتیکی *با درنظرگرفتن 
داده‌های ۱۰ سال گذشته در دو حالت فصلی و ماهیانه استفاده شده و 
طبق معیارهای سنجش و ارزیابی مدل‌سازی» مدل بهتر انتخاب 
گردید. سپس برای هرسال زراعی ۷۱۰۰ سناریو برداشت مختلف با 


باتوجه‌به میزان تغذیه آب‌های زیرزمینی ناسی از بارندگی و میزان 
برداشت در سال‌های گذشته تعیین شده است. به‌منظور حداکثرسازی 
درآمد کشاورزان متناسب با هر سناریو برداشت نیز از برنامه‌ریزی خطی 
استفاده شد. در نهایت عملکرد سه نوع تابع جریمه ساکن, پویا و پوبای 
طبقه‌بندی شده در دو الگوریتم تکامل تفاضلی و ازدحام ذرات برای 
رسیدن به جواب‌های به‌دست‌آمده از برنامه‌ریزی خطی مورد بررسی 
قرار گرفت. شکل ۲ روند کلی پژوهش را نشان می‌دهد. 

به‌منظور مدل‌سازی سری زمانی بارش با استفاده از الگوی 
۸ از نرم‌افزار 10 ۳۷16۷5 و برای مدل‌سازی با روش 
برنامه‌ریزی ژنتیکی از نرم‌افزار 4 100626۳0۲00015) استفاده شسد. 
همچنین برنامه‌ریزی خطی و الگوریتم‌های تکامل تفاضلی و ازدحام 
ذرات به‌صورت برنامه‌نویسی در نرم‌افزار 120185 قتفضا1 1۷۸ 
در دشت سیلاخور نیز با استفاده از نرم‌افزار ۲۲۲۷۷۸۲ به دست آمد. 

در این پژوهش داده‌های مربوط به بارندگی دشت سیلاخور از 
شب رکت آبمتطفهای اشتان لرستان و اطلاعات مربوظ ره محسو لانت 
کشاورزی از ادارات کش‌اورزی شهرستان‌های بروجرد و دورود. 
آمارنامه‌های کشاورزی منتشر شده توسط وزارت جهاد کشاورزی و 
بانک مرکزی به‌دست‌آمده است. 


(90ظ) ممتاه2تصتانن0 نج 5۳ مامتات۲2 -6 
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۳۵۲ نشریه آب و خاک. جلد ۲۵. شماره ۲. مرداد - شهر یور ۱۴۰۰ 


۳ 


نت 


۳ 


شکل ۱- موقعیت جغرافیایی دشت سیلاخور در استان لرستان و ایران 
۵ 200 ۵۲0۲1۳۱6۵ ماعع ۱۵ رز نها سمجاماله ۵۶ ممتایم۱ آهنومج۵و6۵) -1 ۲12۱۳۲۵ 


يارندگي ۱۰ سا اخیر جریمه ساکن 
اعدا با آه ااعکصس۳ ۸۸ ۱ 
۱۱ 


تعبین الگوی کشت با استفاده از برنامه‌رپزی خعلی امقیدسازی مسله جریمد پوپا 
۳ و اس داوم ۱ وتداصای ۵( 0 ۷ 9 


مدل 5۸11۳۸ 
54۱۱۷۱۸ 


مدل برنامه‌ریزی زنتیکی 


اب۲۱ 4:۳ 


3 الگوریتم ازدحام ذرات جریمه پویای طبقه‌بندی شده 
انتخاپ بر اسلس ]۳45 و 7 ۱ 1۳50 واسی عنصصور<ا ف‌تافعم) 
۷ ۵ ۸۳/۹۲ عم اععدط مونءاع۹ 


پبش‌بینی بارتندگی سه سال اینده 
۵ ۱۳۰۰ اد ما مد اعسی۳] الواص۲ 


# بررسی الگوی کشت بهیته 

٩‏ مقایسه عملکرد الگوریتم‌ها برای رسیدن په پاسخ پرنامه‌ریزی خعلی 
ام متا تلاعت اهتاوه که مدناسیینای تسا 9 

۳ ها ترا کم 10۵ ۱ 000۳۱۲۴ ۰ 


تخمین تقذیه اشی از پارندگی 
ااصا وا عم یه ماع 


تعیین ۱۰۰ سناریو برداشت برای هر سال 
0۵۱00 ۱۵ ۵۲۵۲۸۱۷۱۸( 
کف ۷ تاویت ۲ 


شکل ۲- روندنمایی کلی مدل‌سازی پزوهش 


۵06۵ )م۲۵۵۵ 0۲ 060 [۵۵۵۲۵) -2 ما۲ 


مدل مبانگین متحرک خود هم‌یسته بکپارچه فصلی پیش‌بینی سری زمانی فصلی استفاده می‌شود (۵). مدل ٩۸/07۸‏ در 
(5۸1۸) مقایسه با سایر مدل‌ها به تعداد داده‌های کمتری نیاز داشته و رصدی 


فلا ۵40/۸4 یک مان موم خی آنست که از آن برآن کلی بر روی تمامی داده‌ها دارد؛ همچنین این مدل به دلیل دارابودن 


کماسی و همکاران. بهینه‌سازی الگوی کشت در راستای مدیریت منابع آب‌های زیرزمینی با استفاده از الگوریتم‌های 1 و ۳50 ۳۵۳ 


ماهیت فصلی, از توانایی بالایی در مدل‌سازی پدیده‌های طبیعی مانند 
بارندگی برخوردار است. شکل عمومی مدل 5۸77/۸ به‌صورت 
و(5۸/۲/۱/۸)۵,,(<)۳,۳,9 بوده که بخش دوم آن مربوط به 
سری زمانی فصلی است. رابطه ۱ فرم کلی این مدل را نشان می‌دهد. 
5 5 
6( )6( 8) ,0- ,۱۷ ۳ظ) ,۵( ظ) ,4 


۷, <)1-8 (۳ )1-8" | 7 


1 


(ظ),۵ و (8) ,0 به‌ترتیب چندجمله‌ای‌های خودهمبسته و 


از (*ظ),9 ۰ ۰( *8) 6 
به‌ترتیب چندجمله‌ای‌های خودهمبسته و میانگین متحرک فصلی 
هستند. همچنین 8 عملگر وقفه. و 4 به‌ترتیب مرتبه عمل‌گر 
خودهمبسته و میانگین متحرک غیرفصلیء ۴ و ۵0 به‌ترتیب مرتبه 
عملگر خودهیته و ینکن مرک فصلی» 9و یه تیب مره 
تفاضل گیری غیرفصلی و فصلی و ,» دنباله متغیرهای تصادفی مستقل 
و هم‌توزیع با میانگین صفر و واریانس 0 است. مدل 5۸/0۸/۸ 
مبتنی بر سه مرحله شناسایی, ارزیابی و پیش‌بینی می‌باشد (۴ و ۲۷) 


برنامه‌ریزی ژنتیکی 

برنامه‌ریزی ژنتیکی کاربرد بسیاری در مشسخص کردن ساختار 
مشسخص کردن ساختار اساسی پدیده‌های طبیعی یا مصنوعی به‌منظور 
مدل‌سازی عددی آن‌ها استفاده می‌شود (۲۶). خروجی‌ها در مدل 
می‌سوند که اين امر موجب سپولت بکارگیری آن‌ها در کارهای 
اجرایی می‌شود. به بیان دیگر نتایج این مدل کاملا صریح بوده و به 
سادگی توسط سایرین قابل درک و استفاده می باشد. این روش قادر 
است ساختار مدل و پارامترهای آن را بهینه‌سازی کند. برنامه‌ریزی 
ژنتیکی به‌صورت یک تکنیک برنامه‌ریزی خودکار عمل می‌کند؛ 
به‌طوری که جمعيتي برای ایحاد فرزندانی اصلاح شده و به‌منظور 
ژنتیکی به‌جای سلسله ارقام دودوئی» برخلاف الگوریتم ژنتیک روی 
ساختار درختی فرمول‌ها عمل می‌کند. مجموعه توابع و ترمینالها 
ساختارهای درختی را تشکیل می‌دهند (۱۳). 


تغذیه آب‌های زیرزمینی ناشی از بارندگی 

برای محاسبه تغذیه آب‌های زیرزمینی ناشی از بارندگی پیش‌بینی 
شدهء میانگین وزنی بارندگی برای هرسال از طریق رابطه ۲ محاسبه 
شده و با استفاده از درصد تغذیه از بارندگی در طولانی‌مدت. تغذیه سه 
سال آینده تخمین زده شد. روابط ۲ تا ۶ روند کلی محاسبات را نشان 


می‌دهند. 


مور رشن 


7 ۳" 

۳ 2 كثِ_ 
11 1 

57 2-۸ ۴ 

۷۸ - ۲-7 )۵( 

(۶) سک 10 بو 
1 


که در آن‌ها 1 میانگین موزون بارندگی منطقه در سال ام 
(میلی‌متر). رن و مرک به ترتیب بارندگی سالیانه ارتفاعات و دشت 
در سال ام (میلی‌متر. ,4 و م۸ به ترتیب سطح ارتفاعات و 


دست (کیلومتر مربع)» 1 و 7 به ترتیب میانگین طولانی مدت 
بارندگی سالیانه دشت (میلی‌متر) و تعداد سال‌های مورد بررسی» تل. 

7 و ۸ به‌ترتیب میانگین تخلیه» افت و تغذیه طولانی مدت 
(میلی‌متر) 4 و » به ترتیب افت متوسط سطح آب (متر) و ضریب 
ذخیره. ۶ سطح تیسن (کیلومتر مربع) و ۵ درصد تغذیه از بارندگی 
در مقیاس طولانی مدت می‌باشد. در نهایت ارتفاع و حجم تغذیه برای 
هرسال ۶+1 طبق روابط ۷ و ۸ محاسبه می‌شود. 

)۳ زب < ,]1 

۸ رش ی 10*3 پر 

که +: 1 میانگین وزنی بارندگی سالانه دشست در سال ۶+1 ام 

بر حسب میلی‌متر و مر 11 و مر نیز به‌ترتیب ارتفاع و حجم تغذیه 
بر حسب میلی‌متر و متر مکعب می‌باشند (۱۷). 


برنامه‌ریزی خطی 

برنامه‌ریزی خطی پا همان بهینه‌سازی خطی از روش‌های ریاضی 
انش کب فده از آن رات ففتاز پیش یه با که رک تنم 
خطی با تعدادی قید به‌صورت نامعادله راء روی یک چندضلعی محدب 
پیدا کرد. در واقع قیدها به‌صورت نامعادله همان نمایش نموداری 
حتتضامی ,شحو هید ازتباظ آسازم با کاون ک ماخ ارهای اف 
و گستردگی کاربرد این روش در علوم مختلف از جمله مزایای 
برنامه‌ریزی خطی بشمار می‌روند. در این پژوهش برای هر سناریو 
برداشت از آب‌های زیرزمینی با استفاده از برنامه‌ریزی خطی حداکثر 
درآمد کشاورزان تعیین شد. در این راستا تابع هدف مطابق رابطه ٩‏ در 
نظر گرفته شد. 

1 
0 ر۸(ر- رظر 0۲ 2 < 0۵6۷ 
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۴ نشریه آب و خاک. جلد ۰۳۵ شماره ۳. مرداد - شهریور ۱۴۰۰ 


که دراین‌رابطه ۸۷ سود خالص کشاورزان در سناریو مورد بررسی 
(ریال» 17 عملکرد محصول ز ام (کیلوگرم بر هکتار؛ ز] و ر) به 
ترتیب قیمت هر کیلوگرم محصولژ ام و هزینه تولید هر هکتار 
محصول ام (ریال4 ۸ سطح زیر کشت محصول رام (هکتر) و # 
تعداد محصولات می‌باشد. 

همان طور که در رابطه ۱۰ نشان‌داده‌شده است. آب و زمین در 
دسترس به‌عنوان مهم‌ترین قیدهای مدل تعریف می‌شوند. همچنین 
سطح زیر کشت هر محصول به حداقل ۲۰/۶ و حداکثر :۱۷۰ سطح 


زیر کشت فعلی محدود شده است. 
> ۸ 2۲ 
ا< ز 
1 
۸۲ > ر۸ر 
[<ز 
.> رشک 034 


دراین‌رابطه 1۷ و ۷ به‌ترتیب آب موردنیاز محصول ام و کل 
آب در دسترس در هر سناریو برداشت (متر مکعب)» ز۸ و ۸۶ به‌ترتیب 
سطح کشت بهینه محصول [ام وکل سطح قابل کشت دشت سیلاخور 
(هکتار) و ,۸ سطح زیر کشت فعلی محصول ام (هکتار) است. 
الگوریتم تکامل تفاضلی 

یک الگوریتم بهینه‌سازی عمومی علاوه بر داشتن قابلیت استفاده 
آسان باید به‌گونه‌ای طراحی گردد که بطور قابل اطمینان» به مقدار 
بهینه واقعی همگرا شود؛ همچنین زمانی که رایانه برای جستجوی 
راه‌حل صرف می‌کند نباید زیاد باشد. بنابراین یک روش عمومی و 
کارآمد بهینه سازی باید کاملا ساده دارای استفاده آسان» قابل اعتماد و 
سریع باشد. الگوریتم تکامل تفاضلی دارای چنین ویژگی‌هایی است که 
از زمان ایجاد آن در سال ۱۹۹۵ به‌عنوان یک بهینه‌ساز پرکاربرد مورد 
استقبال قرار گرفت (۲۲). اين روش برای شروع کار تنها به تنظیم سه 
پارامتر ۷۳ اندازه جمعیت" ۲ وزن جهش" و 2) احتمال انحام 
بازترکیبی يا تقاطع" نیاز دارد. الگوریتم تکامل تفاضلی از چهار قسمت 
اصلی ارزش‌دهی» جهش,» تقاطع و انتخاب تشکیل شده است .)۲٩(‏ 


الگوریتم ازدحام ذرات 

الگوریتم ازدحام ذرات از رفتار گروهی پرندگان و ماهی‌ها الهام 
گرفته شده است. پرندگان و ماهی‌ها با استفاده از تحربه موفقیت‌آمیز 
خود و اعضاء همسایه, به صورتی که با یکدیگر برخورد نکننده به 
دنبال غذا می‌گردند (۱۱). اين الگوریتم بر اساس همین تقسیم 


ارهز ری ۱۳۹۱۱۱۱۱۱۱۵۵ 
جمنله]ت۱۷۲ -2 


اطلاعات بین اعضاء ایجاد شده و این اصل آن را به یکی از ساده‌ترین 
و درعین‌حال پرکاربردترین الگوریتم‌های فرا کاوشی تبدیل کرده است. 
الگوریتم ازدحام ذرات از دیدگاه محاسباتی سرعت بالایی داشته و به 
حافظه کمی نیاز دارد. 


تابع جریمه 
دو نمونه از مهم‌ترین روش‌های بهینهسازی تصادفی در حل 
مسئله‌های مقید عبارت‌اند از: روش‌هایی که بر اساس جستجوی ناحیه 


قابل‌قبول طراحی شده‌اند 9 روش‌هایی که بر اساس استفاده از تابع 
جریمه مستله را حل می‌کنند (). روش‌های جریمه با اسستفاده از 


عملگرهای جمح شونده و ضرب شونده تابع هدف جدیدی ایجاد کرده 


شونده و ضرب شونده را نشان می‌دهند. 


۲ 1 اک 
(۱۱ 1 ِ ِ ۱ 5( ات 
یوت ۱۱۱ ۱ 5( // ۳ 5( 18 
0 ۲ ]1 ) 1 0 7 


تا ۱۱۱۱9 ۳۹ 7 «( ) / 

( )و عملگر جریمه و 7 ناحیه پاسخ‌های امکان‌پذیر است که 

درصورت خارج‌شدن پاسخ‌ها از این ناحیه. تابع هدف به مقدار (<)ه 

جریمه خواهد شد (۲۲). در مسئله‌های کمینه‌سازی در صورت تخطی 

از محدودیت‌های تعیین شده در حالت جمع شونده ( )2 مقداری 

مثبت و در حالت ضرب شونده مقداری بزرگ‌تر از یک به خود می‌گیرد. 

طبق طبقه‌بندی دیگری» توابع جریمه به دودسته اصلی ساکن و 

پویا تقسیم می‌شوند. توابع جریمه ساکن از مقداری ثابت در طول 

فرایند بهینه‌سازی استفاده می‌کنند در مقابل در توابع جریمه پویاء 
جریمه‌ها در طول فرایند اصلاح می‌شوند. 


جریمه‌های ساکن 

در این نوع جریمه» مقادیر جریمه به شماره نسل‌ها وابستگی ندارد 
و برای پاسخ‌های غیرقابل‌قبول یک جریمه ثابت در نظر گرفته 
می‌شود. هماییفر و همکاران (۸) نوعی تابع جریمه ساکن را ارائه دادند. 
در این روش پس از تعریف .7 گروه تخطی برای هر قید. ضریب‌های 
جریمه نیز توسط کاربر تعریف شده و سپس جمعیت تصادفی با استفاده 
از پاسخ‌های قابل‌قبول و غیرقابل‌قبول تولید می‌شسود. هر یک از 
پاسخ‌ها با استفاده از رابطه ۱۳ ارزیابی می‌شوند. 


3- 01 


کماسی و همکاران. بهینه‌سازی الگوی کشت در راستای مدیریت منابع آب‌های زیرزمینی با استفاده از الگور بتم‌های ۱۳ و ۳0 ۳۵۵ 


0 
[< | 
که 7 به‌عنوان تعداد قیود» رز اهر ضریب جریمه وابسته 
به قید ام و گروه تخطی ام و («اره به‌عنوان قید نامساوی ام 
است. هماییفر و همکاران با استفاده از رابطه 0 > 6-[ ) ره قیود 
نامساوی را به قیود مساوی تبدیل کردند که ۶ یک عدد مثبت کوچک 


است. میخالویچ (۱۶) نشان داد که کیفت پاسخ‌ها حساسیت بسیاری به 
مقادیر این مولفه‌ها دارند. 


جریمه پوت 
جونز و هوک )٩(‏ رابطه پویای ۱۴ را برای ارزیابی پاسخ‌ها ارات 
دادند که مقدار جریمه‌ها با تولید نسل‌ها افزایش می‌یابد. 


1 (۵) 5+( ۶[ مه 


۱ 
۱ 


(م) ,هب2 +(۸۵۶)۲<-(2) ٩0۷‏ 
+< ز <ز 
> ۱ ,0>( )> 0 
(۱۶ تا 9 ۱9/۹۱۱۵ ۱ 0( 8 -[) رم 
( > ز >۱+ و هک[ د) راکو 0 - 
06 ۱ 0 را -(): 


هستند (۲۲). 
جریمه پویا طبقه‌بندی شده 

در این نوع جریمه» گروه‌هایی از تخطی تعیین شده و جریمه 
پویای هر پاسخ باتوجه‌به مقدار تخطی از قیود و شماره نسل تعیین 
می‌شود. رابطه ۱۷ صورت کلی این گونه توابع جریمه را نشان می‌دهد. 

(۱۷) 1۳ ک « () ۳( )+( ) (<-() ۲ 

به‌عنوان مقدار جریمه اصلاح‌شده پویا و / به‌عنوان شماره تکرار جاری 
الگوریتم در نظر گرفته می‌شود. (170 نیز به‌عنوان ضریب جریمه 
مطابق رابطه ۱۸ تعریف می‌شود (۳۱). 


(۴(/]5۳) ۲۱۲۵۲ معتمنوو5 مهع]۷( ۵0۲ -1 


۵ ۲( ,و( ) 2۹ 10 


که 111,... ,1 < ۱,۹( 0( ,0ص ‌ 0( ,۵ است و به‌عنوان 
تخصیص چند مرحله‌ای» ((:7)0 توان تابع جریمه و ():8 تابع فیود 
مسئله می‌باشد (۲۱). 


نتایج الگوی کشت بهیبه 

به‌منظور پیش‌بینی بارندگی سه سال آینده. نخست سری زمانی 
بارش ۱۰ سال گذشته با استفاده از الگوهای سری زمانی 5۸۸ و 
برنامه‌ریزی ژنتیکی در دو حالت فصلی و ماهیانه مدل‌سازی شده و با 
استفاده از معیارهای سنجش و ارزیابی ۱/۸/55 و 96" مورد ارزیابی 
و مقایسه قرار گرفتند. نتایج حاصل از این رزیابی که در جدول ۱ درچ 
شده است بیانگر برتری مدل برنامه‌ریزی ژنتیکی در حالت ماهیانه بر 
ایض بط ی باشزه 

پس از انتخاب مدل برتره بارندگی سه سال آینده زراعی پیش‌بینی 
شد. سری زمانی بارش واقعی دشست سیلاخور در ده سال گذشته و 
سری زمانی مدل‌سازی شده با استفاده از برنامه‌ریزی ژنتیکی در شکل 
۳ نشان داده شده است. همچنین میزان بارندگی پیش‌بینی شده برای 
سه سال آینده نیز در این شکل مشخص شده است. 

در گام بعد تغذیه آب‌های زیرزمینی ناشی از بارندگی پیش‌بینی 
شده در دشت سیلاخور برای سه سال زراعی آینده تخمین زده شد. 
همان‌طور که در جدول ۲ نشان داده شده است. روند بارش و به دنبال 
آن میزان تغذیه در سه سال آینده به‌صورت کاهشی می‌باشد. 

افت تراز آب‌های زیرزمینی در دشست سیلاخور, لزوم تغذیه آن در 
سال‌های آینده را می‌رساند. از طرفی حدود /۵۰ برداشت از آب‌های 
زیرزمینی در این دشت صرف کشت محصولات باغی» مصارف صنعتی 
9 آشسامید تن شده و /۵۰ مابقی صرف کشت محصولات زراعی 
موردمطالعه در این پژوهش می‌شود؛ بنابراین با در نظرگیری حاشیه 
اطمینان» /۴۵ حجم تغذیه در هرسال به‌عنوان حداقل برداشت تعیین 
جو» شلتوک و ذرت علوفه‌ای نیز معادل ۱۴۰/4 برداشت در سال ۰۱۳۹۴ 
مختلف برای برداشت از آب‌های زیرزمینی با فواصل مساوی» تحت 


(1۳0) اصمت‌اهم ومتامصتصه)عب[ -2 


۶ نشریه آب و خاک. جلد ۰۳۵ شماره ۲. مرداد - شهریور ۱۴۰۰ 


جدول ۱- ارزیابی و مقایسه مدل‌های 54107114 و 6۳ 
5 0۳ 21۱0 5۸/۲۸/۸ 0۶ جرمعتوجرجصی 0صه صم6ه 1۲21 -1 12016 


+( مدل (۷۲0061) 
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1 ۲ 
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100 
140 
120 
100 


7 
1 
2 
۳ 
۳ 
۳ 


س حی ص ی کی ای وی 


سری زمانی مدل سازی شده 
65 11۳16 ۷۲006160[ 


و و ات ی وی ای 
زمان (ماه) 
(طاحمص) مصطر]" 


بارش (070) 
۱ 
0 
5 


سری زمانی واقعی 
65 01۳26 16۵1 


شکل ۳- سری زمانی بارش واقعی و مدل‌سازی شده با استفاده از برنامه‌ریزی زنتیکی 
کم عصنوه ۱۵06۱60 0و عمنهی مصت (اهصنه: (هع1 -3 مساعز۲ 


جدول ۲- میزان بارندگی پیش‌بینی شده با استفاده از برنامه‌ریزی ژنتیکی و تغذیه آب‌های زیرزمینی ناشی از آن 
۵۵۵۵ 266۳ )ماع 0ج ط مامت ااقلصنهد 0عاعهعع۲0۲6 -2 م1201 


حجم تغذیه ار تفاع تغذیه 
م۷0 6۵۵2۵۵ . )طاعتهط معبمملع1 
0۷ رصم 
1304 1517 
127 153 
125 12.6 


حداکنر درآمد کشاورزان و الگوی کشت منطقه نیز برای هر 
سناریو برداشت با استفاده از برنامه‌ریزی خطی با محدودیت‌های کل 
زمین در دسترس برای کشاورزی و میزان برداشت تعیین شده در هر 
سناریو به دست آمد. همچنین مساحت کشت هر محصول به حداقل 
و حداکثر /۱۷۰ سطح زیر کشت فعلی محدود شد. قیمت هر 


بارندگی پیش‌بینی شده 
سال 
ار ۱ 
291۳ 
رصص) 
1617 1399-00 
455 1400-01 
148 1401-02 


شده و سود حاصل از کشت هر محصول با بررسی سال‌های گذشته 
به‌صورت نسبتی از قیمت محصول در هر سال محاسبه گردید. 
همان‌طور که شکل ۴ نشان‌داده‌شده است. با افزايش برداشت از 
آب‌های زیرزمینی» درآمد کش‌اورزان نیز افزایش می‌یابد؛ اما در 
برداشت‌های بیشتر از ۲۲۳/۵ ۲۲۵/۲ و ۲۲۵/۱ میلیون متر مکعب به 
ترتیب برای سال‌های زراعی ۱۳۹۹-۱۴۰۰ ۱۴۰۰-۱۴۰۱ و ۱۴۰۲- 


کماسی و همکاران. بهینه‌سازی الگوی کشت در راستای مدیریت منابع آب‌های زیرزمینی با استفاده از الگوریتم‌های 1 و ۳50 ۳۵۷ 


۱ محدودیت کل زمین قابل کشت مانع از افزايش سطح زیر کشت 
همچنین منحنی درآمد هرسال از پنج خط تشکیل شده است که در 
طول هر پاره خطء سطح کشت یک محصول افزايش پیدا کرده و 
باعث افزایش درآمد می‌شود. در طول پاره خط پنجم نیز سطح کشت 
تمام محصولات ثابت باقی‌مانده که به دنبال آن درآمد نیز ثابت 
می‌شود. همچنین به دلیل در نظرگیری تورم در تخمین قیمت سه 
سال آینده محصولات کشاورزی؛ منحنی درآمد در طی سال‌های 
متوالی افزایش‌یافته است. 

شکل ۵ مساحت کشت بهینه حاصل از برنامه‌ریزی خطی را در 
سناریوهای برداشت مختلف برای سال ۱۴۰۰ نشان می‌دهد. 
ترتیب تغییر مساحت کشت محصولات بیانگر آن است که اولویت 
افزایش سطح کشت یک محصول همزمان افزایش میزان برداشت» 
رابطه مستقیمی با به‌صرفه بودن آن محصول از نظر سود خالص و 
مصرف آب دارد؛ به بیان دیگر با مصرف مقدار مشخصی آب به ترتیب 
افزایش برداشت آب برای دو سال ۱۳۹۹-۱۴۰۰ و ۱۴۰۱-۱۴۰۲ نیز 
به همین صورت است. 

با مقایسه الگوی کشت بهینه حاصل از بنامه‌ریزی خطی و الگوی 
کشت بپینه فعلی که در شکل ۶ نشان داده شده است می‌توان 
دریافت» در صورت عدم تغییر میزان برداشت از آب‌های زیرزمینی و 


۰. 1401-02 


00-1 
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برداشت آب زیرزمینی (میلیون مترمکعب) 
(۷۲۵) ممتاماتما۱ ماد همین 


پیروی از الگوی کشت بهینه» درآمد کشاورزان منطقه ۱۸/۶ افزایش 
می‌یابد؛ به صورتی که» سطح کشت محصولات گندم و ذرت علوفه‌ای 
یت یه مخ کمت افیلی فا سوه مادم رنف ده آفر یش بافیة 
و سپس به‌منظور حداکثرسازی تابع هدف به دلیل صرفه اقتصادی 
محصول شلتوک نسبت به جو سطح کشت آن افزوده می‌شود. 


نتایج بررسی الگوریتم‌ها و توابع جریمه 

برای سه سال آینده عملکرد دو الگوریتم تکامل تفاضلی و ازدحام 
ذرات برای حداکثرسازی درآمد کشاورزان و رسیدن به جواب‌های 
به‌دست‌آمده از برنامه‌ریزی خطی بررسی شد؛ در این راستا به‌منظور 
نامقیدسازی مسئله تعیین الگوی کشست. محدودیت‌های زمین 
قابل کشت و آب قابل برداشت در هر سناریو به‌صورت توابع جریمه 
ضصرب شونده تعریف شسدند. همچنین تعداد جمعیت و ضرایب 
عملگرهای تقاطع و چهش در هر دو الگوریتم یکسان در نظر گرفته 
سال ۱۳۰۰-۱۴۰۱ نشان‌داده‌شده است. به‌طورکلی الگوریتم ازدحام 
ذرات در این مسئله بسیار سریع‌تر به جواب می‌رسد. به طور میانگین 
تعداد تکرارهای الگوریتم ازدحام ذرات برای حل هر سناریو در این 
متسه ۳۲ تاه فخرارهای الگور کال اخیلی ای لگ 
کال اسان تا تضیر هن انط متفر هقی ازهای تلازم راخ 
رسیدن به جواب تغییرات محسوسی از خود نشان می‌دهد. 
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شکل ۴- درآمد خالص کنساورزان متناسب با برداشت آب 
اهامای مه ۵ حمتا تامحز عم ره ۵ مصرمم‌ه ا۵ 11۳6 -4 مساع۲1 


۳۸۸ نشریه آب و خاک. جلد ۰۳۵ شماره ۳. مرداد - شهریور ۱۴۰۰ 


ذرت علوفه ای ۰ ۳ < شلتوک ۰ 
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برداشت آب زیرزمینی (میلیون مترمکعب) 
(۷۲۱۵۵) ممتاهاتماند۴ ماد صنمین 


«ممه 3 ۷۷۵۵۲ 


مساحت کشت (هزار هکتار) 


(دعتهامع۲۱ 103) معته ومتن 


50 710 90 110 130 


بوک زرح 


شکل ۵- الگوی کست بهینه محصولات کشاورزی در سال ۱۴۰۰-۱۴۰۱ 
2021-72 ۱ ماهس اوه ۵۴ )اه موم اقصتامون0 -5 ۲120۳6 


۵ 00۵ میت 


لح کشت بهی: 
۵ 0۵ افصتاون0 


ذرت علوفه ای شلتوک 


مت 1۹16 


21 


سطح کشت (هزار هکتار) 
ی 0 00۵۵ 
مر 


جو گندم 
7 ۵ 


شکل ۶- الگوی کشت فعلی و بهینه با میزان برداشت آب یکسان 
0 ۷۵۲۵ 0 )امه موی م6 6 ماه عصامممن آفصتانون 0ص )همست -6 م۲120 


منظور از تغییر شرایط مسئله, تغییر محصول انتخابی برای افزایش 
سطح کشت با افزايش برداشت است. پنج ناحیه متناظر در شکل‌های 
۴و ۵ که شرایط مسئله در آن‌ها تغییر می‌کند نیز در شکل ۷ 
نشان‌داده‌شده است. درنظر گرفتن یک تابع جریمه یکتا برای تمام 
مسائل بهینه‌سازی به دلیل حساسیت بالای کیفیت پاسخ‌ها به 
پارامترهای جریمه» غیرممکن است؛ اما طبق نتایج مندرج در شکل ۸۷ 


اعمال توابع جریمه پویای طبقه‌بندی شده پوبا و ساکن به ترتیب 
بیشترین تأثیر را در کاهش تعداد تکرارهای الگوریتم برای رسیدن به 
پاسخ بهینه حاصل از برنامه‌ریزی خطی ایفا می‌کنند. روند کلی تعداد 
تکرارهای دو الگوریتم برای رسیدن به پاسخ‌های به‌دست‌آمده از 
برنامه‌ریزی خطی در دو سال دیگر نیز به همین صورت است. 


کماسی و همکاران. بهینه‌سازی الگوی کشت در راستای مدیریت منابع آب‌های زیرزمینی با استفاده از الگوریتم‌های 1 و ۳50 ۳۵۹ 
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شکل ۷- تعداد تکرار الگوریتم‌های تکامل تفاضلی و ازدحام ذرات با سه نوع تابع جریمه برای سال ۱۴۰۰-۱۴۰۱ 
2021-02 ۲۵۲ عصمتاص پاامصون ان ما م66 ۱۱۵ عصط)تدمعاه 0فظ فصو را ۵ عصمت )1 گ۵ تمجاصصت ۱ -7 ۱۵ع۲۱ 


تتایج مندرج در جدول ۳ حاکی از آن است که الگوریتم تکامل 
تفاضلی همان‌طور که نسبت به تغییر شرایط مسئله واکنش نشان 
می‌داده در اینجا نیز با تغییر تابع جریمه از حالت ساکن به پویای 
طبقه‌بندی‌شدهء تعداد تکرارهای آن برای رسیدن به جواب برنامه‌ریزی 
خطی به طور میانگین ۱۱/۶ کاهش می‌یابد؛ در مقابل الگوریتم ازدحام 
ذرات نسبت به تغییر تابع جریمه واکنش محسوسی نشان نداده و 
تکرارهای آن به طور میانگین تنها ۲۳/۸ کاهش می‌بابد. 

در شکل ۸ تاریخچه زمانی فرایند بهینه‌سازی درآمد خالص 


کشاورزان» در ۱۰۰ تکرار ابتدایی دو الگوریتم. نشان داده شده است. 
این نتایج متعلق به حداکثرسازی درآمد در سال ۱۴۰۰-۱۴۰۱ و 
متناسب با سناریو برداشت ۱۵۰ میلیون مترمکعب می‌باشد. همان طور 
کهتقر تک ی اس گرفه گوزش امن سناضتان با اعماز. 
تابع جریمه ساکن پس از ۲۶ تکرار و با اعمال تابع جریمه پوبای 
طبقه‌بندی شده پس از ۱۷ تکرار به پاسخ بهینه همگرا می‌شود که 
بیانگر تآثیر توابع جریمه در فرایند بهینه‌سازی این الگوریتم می‌باشد. 
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جدول ۳- میانگین تعداد تکرار الگوریتم تکامل تفاضلی و ازدحام ذرات برای تعیین الگوی کشت بهینه در ۱۰۰ سناریو 
۹ 100 1۱ )وج عصتم‌ومی اعصتامرن عمط ماع 6۵ عصصطتدمعاه مق هصرع بر اه مصمتامهاز ۵ بحاص عهه۸۲ -3 م1201 
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میانگین درصد کاهش تعداد تکرار الگوریتم نسبت به تعداد 
تکرار با جریمه ساکن 
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شکل ۸- تاریخچه زمانی فرایند تعیین حداکثر درآمد کشاورزان 
۵۵۵۵ 1 م160۱ حصاحص هه معط هصتصتصماع0 ۵۶ وععیم معط ۵۶ رماعنطا م1 -6 مساعز۲ 


الگوریتم ازدحام ذرات نیز با اعمال تابع جریمه ساکن پس از ۱۳ 
تکرار و با اعمال تابع جریمه پویای طبقه‌بندی شده پس از ۷ تکرار به 
پاسخ بهینه میل می‌کند که اين موضوع بیانگر تأثیر کم تغییر توابع 
جریمه در اين الگوربتم می‌باشد. روند کلی تاریخچه زمانی فرایند 
حداکثرسازی درآمد کشاورزان در سایر سناریوهای برداشت نیز برای 
هر دو الگوریتم به همین صورت است. 


نتیجه گیری 
در این پژوهش ۱۰۰ سناریو برداشت مختلف از آب‌های زیرزمینی 
برای سه سال زراعی آینده باتوجه‌به تغذیه ناشی از بارندگی پیش‌بینی 


شده و میزان برداشت در سال‌های گذشته تعیین شد. سپس الگوی 
کشت بهینه با تابع هدف حداکثرس ازی درآمد کشاورزان و 
محدودیت‌های آب 9 رمین در دسترس برای هر سناریو برداشت با 
استفاده از برنامه‌ریزی خطی به‌دست‌آمده و عملکرد دو الگوریتم تکامل 
طبقه‌بندی‌شده برای رسیدن به جواب‌های به‌دست‌آمده از برنامه‌ریزی 
خطی مورد بررسی قرار گرفت. نتایج نشان می‌دهند الگوی کشت 
منطقه بهینه نبوده و با بهینه‌سازی آن, می‌توان درامد کشاورزان را تا 
افزایش داد. همچنین با افزايش برداشت آب. الگوی کشت 


کماسی و همکاران. بهینه‌سازی الگوی کشت در راستای مدیریت منابع آب‌های زیرزمینی با استفاده از الگور بتم‌های ۱۳ و ۳0 ۳۶۱ 


و توابع جریمه نیز حاکی از آن است که در حل این مسئله الگوریتم 
ازدحام ذرات نسبت به الگوریتم تکامل تفاضلی در ۸۴/ سناریوها با 
تعداد تکرار کمتری به پاسخ بهینه می‌رسد. اعمال توابع جریمه 
ظقپری شنده نی یاه تکرارهای الکتری تکام تفای را 3 ۸ 


در نهایت پیشنهاد می‌شود سیاست‌های اقتصادی حاکم بر بخش 
کشاورزی با درنظرگرفتن میزان سود خالص محصولات کشاورزی 
نسبت به میزان آب مصرفی آن‌ها اتخاذ شود که در صورت افزایش 
برداشت آب از منابع زبرزمینی» الگوی کشت به سمت محصولات 


تعداد تکرارهای الگوریتم ازدحام ذرات را //۲ نسبت به استفاده از آب‌بر نرود. 
جریمه ساکن کاهش می‌دهد. 


منابع 
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مصلوا تلهم عصلم۳0۵ عص2تصصتان(0 :2015 .ن1 وملاه مه ر.ق سعلد ره .۱۷.۲۲ مصمصصیان .1 اعدا رک تاصمطاتاظ 
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1245-7۰ :(۱91)11 عطزععمتممظ۴ مه وعتصمطمع۱۷ 0عنام۸ ظا ما۱۵ عمجم .2 فطا ۵۶ ماهاه مطا که بمتناو ۵ 
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مصازممصامصرظ 1ماصمصصمم۱1 ۱ مرج زر ۵۶ امصتامز .ملق متصرنا :رتاو مقمت متافمه همم ۴۱۸۲۵۵ 
امه طفتاوصط از مقتو۳۵ 19) :75-86 :(45)3 و1۵02 0۳ 

۵ 2۷۵۲ طا مادم فعصتونمته فطا ۵۶ ممتاه‌متصصتامون0 .2017 1,۸۲ طمعصقکا مضه ری ت«امطلاط ریگ مهو 
36)4(:557-۰ رلقصتتاهز عطاهعصلم ۴۸ 
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(۳۵۲۹120 1) ,مهم مد معصمععتصمن مفتامملهو معصعااقطه 
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نشریه آب و خاک (علوم و صنایع کشاورزی) 8 60 0ص ۲۷26۵ ۵۶ تمه[ 
جلد ۰۳۵ شماره ۳. مرداد - شهریور ۴۰۰ ص. ۳۴۹-۳۶۴ بو زسه 49-4 .0 ,2021 .۸820-.01 و3 ۷0۰ و35 ,۲۷۵0۱ 


۵۹۵ 6۶۲۰0۱۵۱۱۱۲26۳ ۳0۵۵ )و۲ منونمم 0 صماه2صن6 
م۸۱۵0 ۳9۵0 20 انا مصزو ۱۷۲2۱۵220۵006۵0۵ 


1۷۲, ۵۳29۲ - ۸ 7 
16061۷60: 23-02-01 
66060: 04-05-01 


0۵ 0۶ 1201 0ج ۱۵605 مه مصتفهممصا طازه عجماه رومام عتووعممناو ۵۲ ومهمتتناهعن 18 خمملام00 1۳۸0۲۵ 
متصمجمعع همه لمامممصممت هه فیامتتهه ‏ مق قطا ملفتته. ماه 2 میاه مق فممتامومع ما۳۷ 01 اصعطصمع۳۸۵2۵۵ 
6 2180 مه دوم م۷2 فهت۱0۳۷2 مهم عصتوممین فط ما معصفطه مطا رهام فط مه 200100 ق .عععصمنتمعقهمم 
مطا ۵۵ )م1 امعنل ه فقط مت ,زاتمم فطع ما وتقع۷ ۲۵۵۵۵۲ صد مممممممصمای لهمتاتن اف مطا دتفتت 728۲ عطا 
وم مصا ۶ مم و1 ماه مملوون ما امه ,۲۳۵۵10۲ ما۷۵ 0۴ تممصتافومن اوعع127 مصا وه ۵0۵۲و امتنآتام11ع2 
ما محتصعاعل ما ,لاه فنص ما .قمعه) تمه مج طا عمتجم هه وعمتناموع۲ تماه/ عظع0202 صد م20 ]1۳0001 
تقعصن! منود فعطمه0؟200 ۷۵ عصلفتا و۷62۲ ۲۵۵ظ) ۱670۲ عط ما صتهام عمحصلم11ه صا دمم تمزمصط گم مادم عصلمممنن لقحصتامه 
۰ 15)1 ۱۷۵۲۵-۲۲۵۱۲ 4مه عمتصصصهه۳۳۵ 


۶ صنح معط ۳ مادم عصتمومن لقصتاوه مطا مصتصعاعل ما تعل صا رمرماه افیا فطع صا ءفل‌۱۷]۵۵۵ 0ص2 واهتعم)۱۷2 


الصه ۶ )مجح عمط رع2۷211201 1200 مه تمه۳7 10160 فص ۵0 و۷۵2 ععط) )»رم مط) صا دممطممهصا عتمطنه عصمتصطل تم 
40 تمه فطع مصتصمامل مه لهمنصا ممتاه)امآمونه تمله۳ معط مصتصصعاعل ما ممتماتته 2 وه ۱۷۵6۵0 15 ۲۵۵2۲826 
0 24 واع0مص وعتیعه مصصتا ۸۲۸و بللفکصنه: بفدهع0۲؟ ما تع0ع۵ ص1 تقع۷ طمقه صمته)امانه حصیاصصتته۳ 
0 روع0مص باط‌جمصه مه لقصمقدعه طامه ص وتوع 10 )فق1 فا ۵۶ 0202 فمط) عصتیملنفوم 64و م۷۵ مصتصصهع۲۵ظ 
ممتهل0010 100 رتهع۷ ۵۲00۵ طمهه 10۲ رحقط 1 .60)مع۵1و ۷/۵8 1ع0مظ تمتاهط 2 بدتعاتت .۲۱۱ مه ۳۸9 ما عصتل001ع2 
۶ )صیامصصح مه 0ص 1ص ها 0عمناه مفتهطمع۲ ماه نامع 0 اصیامصصع معط ما عصتل01عع2 60صتصعاع0 ۷/۵۲۵ 866927105 
۰ 01۲6۷10۱۷5 طا ممتاقا نارهم 

۶ نج مط) 1 مهم عمتوومت اقصصتاوه مطا ماع م6 ۱۵۵0 ۷2۵ مصتمه ۳۳0۵ فصن روماد 960000 مط 1 
۶ 01160 1۳6 ,1220 ماصاجته 20 متتقجهعز طمجه ما تمه ماصهاتماه ۵۶ مصمتاقاتصصا مج وعمصممهصا فتمصرنع عصام2تصل(۳۳2 
۳0 ]۳01 مطا 20 ,206 .1.6 ,اصنامی عمط گه صمتاحاگصا حیم)-عصم1 عههعع2۷ معط م) عصتل 20601 0۳0[6۵)60 15 )ما۵0 مهم 
مهن ها کهع موجه صا امنالمص مط ۵۶ ممتعج مطا ۵۶ ممتازمو۵۲0 ه وچ )12اه قه ]ع0مم طمقع ۶ ممتاه تاه عطا 
۵۲ ومتصهه دنا 0ع1اوود1ن) مه رمتصمطرن رمتقاگ 0۶ موجه معط ۵۶ ممصمصصمنه فط ,افص قمع منامانميم فط 
0مصنهاداه فالناوع عمط عبمنطمه م۵ 1۳۷690122460 مه 80ظ م4صه متام اهت)ممنعذن1 رفحصطاتنمعاه و مفصا فصمتام‌صا۴ 
8 ,0۲06659 220101 مان مط) عصیل مباله۷ اصمافجمع 2 وبا فصمتاعصیظ کالقصوم متاهاگ .عصتحصطنهع۱۵زظ تحمصننا ۶۲۵ 
۰ 0۶ ۲تامصتاه مص ون مصومع0 220 وقععمعم مطا عصتتل ۵016160 مه ففص مطا رفصمتامصظ القصوم متصصهصرول ظ1 
8 ۲۵5۵0۵86 مه ۵۶ هاقصهن معط 20 رماع مولج مه فصمتاهام۷1 ۵۶ دومع رسالقموم فعتصصهمرل 11160ووهاه مط 1 
۰ 2010عصعع مط) مه فصمتام‌تتاوع1 عمط ۵۶ ممتاهام۷1 ۵۶ اصتامصصح ما مه عصتعمعع2 0مصتهاع 

ومصمعص1 قمصنق] رومتاماتمان تمه صتمرع مصتفممتمص1 ره اقطا مدای فالناوعر م1 خصمتععیهعزظ م0صد و6لنعع1 
ما 101 وتمامصه مذطایت ممتلاتمظر 225.1 280 ر222.2 ,223.5 صقطا معمم ۵ ممتامااماه مطا مر ۲۱۵۱۷۵ جموهمهه1 2190 
فقط 0صه1 ماطادیه ماما مطا ۵۶ صمتاداتصصنا عطا ۲۵50۵6۵۳۷ ,2022-2023 20 ,2021-2022 ,2020-2021 و۷62 مصنم0۵؟ه 
۰ 1۲۵۳۱۵1 وعصطمعص ومصته] رجمتاهاامآمونه مصلوهم مهد رها م2 مصمتاه تن تعصتا 2162 معط ۵۶ وقهعععص1 مط 0۲6۷۵۵۸۵۵0 
مماه۲ تاه تمصرنا قمع اممتهن مطا مه رمع 20 رعع1 ررماتهه بتقمطاه ۵۶ دوم عنام ماه تانه ما تعقت صا روولظر 
۵ ۱ .5۵۱۲۵۵5 صیامع تعصا رم «متصصع 0ع)وع7 بط و1 عمیه/7 0۶ وتمامجط مزمایت ممئالنمه 142 رصتقام تمصلمان 
6 ما وتمصصرتها ۵۶ مصطممصا مص رهمتاماتمامونه ۷28 آممستن فطع هه رلعنبتاو دومن عنام مفط م۶ ممتاهم عصونومه معط 
۰ 0 186۲6286 ۷1۱1 ۲۵۵108 

0 عط) ۶ دعمتاح‌عها1 که میتی عع2۷62 1۳5 .تماو2؟ طمننجر مصرماطاهعم فتطا فعسعصه حصطتع0ع2۱ 0وظ معط رلهتعصهع و[ 
۱۷۵0۵۵۵۵ )1۱۱1۶1۵۲6 معط ۵۶ فصمتاحتعا1 که تعطامصننه مط) ۵۶ 3860 ور مصماماميم ونطا ما متتقصوعو طمقع 501۷۵ مه 212011 
ص .ومنتقصوعو عطا ۶ه 840 طا تعتاعه 0مصصتمتتهم فقط نارمع( 90ظ مطا رصعاهامم فنطا عم ا۵د مد رالهعن .حصطتت0معل2 
متقاه 240 ,متصصهطول رومتصهم رل 160رووهه عط مه فهمماه‌ها ممصط)تعمعاه امه ظد ععصقمموج ]وه فطا رفعمتامصنظ سالقصهم 
۷ ومتصصه ول 11160ووهاه ما متلقاو حصم ممتامصنظ جالقصوم معط عصتعصقطه اظ ,اهتاممموع۲ مقطمتاعصظ سالقصهم 
عصتصصهوع ۳۳0 بهمصنا مط ممتطمج مه صاتزمعاه ممتتاهظ توتامم۲ع]۱۱۶ معط گم معمتاجمهاز ۵۶ عطاصته مطا مصمتاعصی] 
صفقط مطا ما امه نصوله ۲۵۵۵۲ ]۵ 10 صطتعمع(2 0وظ مطا رافمتاومی ص1 م11 ۵۶ م2۷۵2 صه زا ۲۵00660 و1 ممتانا01و 
۰ ۵01۷ 0۶ 2۷6۲۵۵ مه 0۲ 060۲62860 فصمتاتا۲6۵6 فا فصه مصمتامصنظ پرالقصوم عطا طر 


نع زا۱010 طمااماه ره روه‌تامند متاییمت بواظ مه مصتهعصتمم تماج ۷۷ ۵۶۴ ۲۶۵1۵950۲ )9002و 24 اصمن0 ,1۷۲9۹۵ -2 20 1 
۰ 121 ,له زنا010ظ تفه نون 

(2010.26.11 ۵ تومعصم! تاتمعصظ ممطنه م مت طمصدمتزمن -۶) 

10۵]: 5 


۴ نشربه آب و خاک. جلد ۰۳۵ شماره ۳. مرداد - شهریور ۱۴۰۰ 


۶ 1062866 عط) از مج ,لقتامه )0ظ وا ممتعع۲ مطا ۵۶ مادم عصتمومه مطا تقطا «مدو ناف م1 خممتعسمومن) 

ممتاه تایه عمط متام رها رمفله تلم ماج که ممتاهب نت مط) قهمهوه ومصصر له 1 رطمتاهااماه ما۳2 
و موه متام همتاممع نا گه ممتاممتصیدظ .60و2عععصا هصق ومحرممصاً. فعمصیه رطمنعع۲ عطا گم )22 
ص متامصنط الموج وعتصحطول 1)160ووهاه مطا عصتوتا عقطا ب«مطی موق فصمتاعصیظ رالممهم گم دعمها معط ظ16 وحصطا 212011 


۰ 2000 2۷6ظ صقه حصطاتت0ع21 90ظ فط 


صتقام تمصلمانه رجمتهمتصتام۵ رفصطتمماه متافتسهمطمام۷ رمصتصمصهع ۲۶۵ تم رماه امین :1۵۲۱۷۵۲9 


